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The Bohr hamiltonian, also called collective hamiltonian, is one of the cornerstone of nuclear 
physics and a wealth of solutions (analytic or approximated) of the associated eigenvalue equation 
have been proposed over more than half a century (confining ourselves to the quadrupole degree of 
freedom). Each particular solution is associated with a peculiar form for the V(/3,7) potential. The 
large number and the different details of the mathematical derivation of these solutions, as well as 
their increased and renewed importance for nuclear structure and spectroscopy, demand a thorough 
discussion. It is the aim of the present monograph to present in detail all the known solutions 
in 7— unstable and 7— stable cases, in a taxonomic and didactical way. In pursuing this task we 
especially stressed the mathematical side leaving the discussion of the physics to already published 
comprehensive material. 

The paper contains also a new approximate solution for the linear potential, and a new solution for 
prolate and oblate soft axial rotors, as well as some new formulae and comments, and an appendix on 
the analysis of a few interesting numerical sequences appearing in this context. The quasi-dynamical 
SO(2) symmetry is proposed in connection with the labeling of bands in triaxial nuclei. 
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I. INTRODUCTION 

More than half a century has elapsed from the publication of one of the milestones of nuclear physics "The coupling 
of nuclear surface oscillations to the motion of individual nucleons" [1] on the danish journal "Matematisk-fysiske 
Meddelelser" by Aage Bohr in 1952. The cited work contains the foundations of the collective model, developed fully 
in a second important work in collaboration with Ben R. Mottelson [2] , and especially contains the first derivation of 
the subject of this monograph: the Bohr hamiltonian, also called collective hamiltonian or , sometimes, Bohr-Mottelson 
hamiltonian(i7s). This operator contains kinetic as well as restoring potential terms in some set of variables describing 
the extent of deformation of the nuclear surface. The present paper, far from being a comprehensive review (that are 
numerous in the specialized literature [3, 4]) of the many physical topics, successes and applications of the collective 
model, is rather meant to provide the reader with an enumeration and mathematical discussion of the analytic and 
approximate solutions of the stationary eigenvalue equation = with various form of the potentials and 

referring to different physical situations. Other developments originated from the collective model (as the Frankfurt 
model, the Baranger-Kumar model or the symplectic model) will not be discussed here. 

The need for such a taxonomic compendium is, to our view, justified by the recent proliferation of new solutions 
that makes this topics exciting and, at the same time, significantly extended. The new solutions have received a 
considerable attention, both from theoreticians and from experimentalists, and application of these ideas or survey of 
spectroscopic data are constituting a very important research line. It is also our belief that these studies have reached 
a considerable degree of completeness and that a classification of the subject is therefore timely To corroborate our 
views upon the renewed role of the collective model as a driving force for new research, we mention that an interesting 
paper, that contains some sections dedicated to the Bohr hamiltonian, has recently been published [5] (See section 
VF). 

The problem of the solution of the Bohr hamiltonian has been mainly tackled in two ways: direct solution of 
the second order differential equation and use of algebraic techniques to exploit dynamical symmetries. The latter 
approach has also been used as a source of 'labels' for various solutions of distinctive importance. A novel interest 
has been raised by the possibility to give analytical solution at the critical point of a shape phase transitions between 
various types of quadrupolc deformed nuclear surfaces [6, 7]. Another method that is worth mentioning by virtue of 
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the insight that one may get is the numerical diagonalization of the problem for complicated potentials that cannot 
be treated analytically. We will also remark on some of these cases. 

After a brief discussion on the historical and scientific foundations of the collective model (Section II), we will be 
concerned with an analysis of the known solutions of the eigenvalue equation for the Bohr hamiltonian for 7— unstable 
(Section III), axial 7— stable (Section IV) and triaxial 7— stable (Section V) cases, with direct solution of the differential 
equation, which we will lay before the reader following, as far as possible, the course of time. Some comments on 
group theoretical techniques and solutions are also given (Section VI). Section VII contains a brief description of a 
few important recent advances that are connected with the main theme of this compendium. 

Writing in a single paper a summary of works and efforts that have been devoted to the study of collective models 
and paying the tribute to every researcher involved in the matter would be a formidable task, and we have decided to 
confine ourselves to the work plan cited above, apologizing for any unintentional omission. We decided to label the 
solutions on the basis of their proposers, mainly to avoid confusion with the corresponding problems in the solution 
of the Schrodinger equation in the configuration space. Where the authors were many we used alternative titles. 

This paper is not only a review, but contains some novelties as the treatment of the linear potential in Section 
III.H, that to our knowledge has never been treated in this context, or, for instance, formula (40). Moreover we 
propose a new solution for prolate and oblate soft axial rotors (Section IV E) obtaining interesting patterns for 7 and 
/? excitations. 

A critical discussion of the various cases is supplemented by few comments that may share light on specific points (as 
for instance in section III 1 2). In the appendix we discuss from a purely number theoretical perspective, the problem 
of repetitions of quantum numbers, that, as far as we know, has never been discussed elsewhere: we give two new 
sequences and we evidence some curious connections with other fields. 



Niels Bohr, father of Aage Bohr and of the atomic theory, is among the proposers of the theory of nuclear surface 
oscillations [9]. The idea was to treat the nucleus as a liquid droplet, and that the fundamental collective modes of 
motion of the surface were linked to nuclear excitations. Let us quote Aage Bohr's words [1] about the liquid drop 
theory: 

According to the liquid drop model, the fundamental modes of nuclear excitation correspond to collec- 
tive types of motion, such as surface oscillations and elastic vibrations. Even if it has not been possible, 
with certainty, to associate observed nuclear levels with particular modes of oscillation, the model gives 
an immediate explanation of the rapid increase of level density with increasing excitation of the nucleus. 

The fifty years elapsed since the time of the exposure of his ideas have contradicted the last sentence: it is indeed 
possible to associate solutions of the Bohr hamiltonian with nuclear spectra, giving an immediate explanation of many 
nuclear properties. The collective model, that was subsequently developed in collaboration with Ben R. Mottelson 
and that has taken the name of Bohr-Mottelson model, should be regarded, as Bohr was warning, as a complementary 
approach to the shell model and should shed light only on some aspects of the still dim complexity of nuclear spectra. 
The innovative idea underlying the already cited work was 

...to consider various properties of a nucleus described in terms of a deformable surface coupled to the 
motion of individual nucleons. This combined model may be referred to as the quasi-molecular model. 

It is strange that the name proposed by the author has been lost. He found that his model bore a resemblance to the 
theory of molecules in the very same way in which the single particle shell model bore a resemblance to the atomic 
theory. 

Bohr's paper was confined to the treatment of a single nucleon interacting with the nuclear surface. The radius of 
the latter was expressed, in polar coordinates, as an expansion in spherical harmonics: 



where Rq is the radius of the nucleus when it has the spherical equilibrium shape. The expansion parameters are 
the coordinates that define a multidimensional space, whose points represent a deformed surface. The requirement of 
reality for the nuclear radius implies 



II. FOUNDATIONS OF THE COLLECTIVE MODEL 




(1) 



(2) 
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The consistency of expansion (1) with the 'atomicity' of the nucleus, i.e. the fact that it is composed by smaller 
particles, the nucleons, is assured by realizing that the idea of a continuous surface loses its meaning if one takes into 
consideration surface elements whose dimensions are comparable with the average distances between the constituents. 
This happens whenever A is larger or equal to the linear dimensions of the nucleus, A 1 / 3 . 

In the following we will confine to the quadrupole degree of freedom (A = 2) leaving aside the discussions about other 
multipolarities: in fact quadrupole states are the most fundamental collective type of low-lying excitations in nuclei 
and deserve a thorough analysis. If only small values of a's are taken into account the quadrupole deformed surface 
is an ellipsoid randomly oriented in space. The five coordinates {a^,^} may thus be mapped onto a set of other five 
variables {cio, 02, #i, #2, $3}, three of which describe angular orientations of the ellipsoid, the other two parameters 
being related to the extent of the deformation. This mapping from the frame of reference R to a second frame of 
reference R! , called 'intrinsic', may be achieved via a unitary transformation: 

2 

a% >v = 2J a 2.^^. y {9i) (3) 

iU=-2 
2 

<*2 t(l = Yl a ^KA i)> ( 4 ) 

ju=-2 

where T>^^ v {9i) are the Wigner rotation functions for spherical harmonics with A = 2. The set of three angles 
describes the relative orientation of the two frames of reference. From now on we will drop the index 2 since, as 
already said, we will deal only with quadrupole deformations. The transformation (4) may be conveniently chosen so 
that the axes of the frame of reference R! coincide with the principal axes of the ellipsoid. In this case 

ai = o_i = and a 2 = a_2 • (5) 

The transformation of coordinates described above is however not unambiguous because of the possible different 
order of labeling of the axes and of the choice of the direction of the versors. This amounts to 48 sets of tranformcd 
coordinates which correspond to the initial set (24 as long as right handed coordinate systems only are considered). 
The symmetry properties of the wave functions will be affected by the arbitrariness of the choice of the transformed 
coordinate system and this fact implies certain symmetry relations. Another substitution is required to reach a new 
set of coordinates {/?, 7,^} called the Hill- Wheeler coordinates, in which the subset {(3, 7} defines a two-dimensional 
polar coordinate system within the five-dimensional quadrupole deformation space, namely: 

a = (3 cos 7 (6) 
a 2 = (/3/ V2) sin 7, (7) 

with ^ I 1 2 = p 2 . The total deformation of the nucleus is measured by (3 in such a way that (3 = represents a 
sphere and [3 ^ is an ellipsoid. The larger the value of (3, the more deformed the surface. The parameter 7 describes 
the deviations from rotational symmetry: whenever 7 = n^, with n £ Z, two of the three semi-axes are of equal 
length. This is depicted in fig. 1. 



To each (solid-line) radius in the polar plot of this figure corresponds an axially symmetric ellipsoid with a well- 
specified symmetry axis and a given character for what concerns prolateness or oblateness. Every point in the areas 
within the radii corresponds to triaxial shapes, that are characterized by three different semi-axes. 

A convenient way to deal with quadrupole surfaces is to rewrite eq. (1) in terms of the new variables. The three 
radii defining the ellipsoid in the body-fixed frame are thus 



Rk = Rc 



5 2 
1 + — (3 cos (7 — -irk) 
47T w 3 ' 



(8) 



with k= 1,2,3 . 

Now we are ready to introduce the Bohr hamiltonian, that is the hamiltonian built with generalized coordinates 
and momenta in quadrupole deformation space: 



n 



(9) 
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FIG. 1: Hill- Wheeler coordinates {/3, 7}. The quadrupole deformed shapes corresponding to j3 = 0.4 and 7 = mr/3 (with 
n — 0, 5) are shown for reference. Different colors are connected with different principal axes of symmetry (green for z, red 
for y and blue for x). 



where B2 (that will be called B m i n the following) and C2 are the mass and stiffness parameters of the liquid drop 
model for the quadrupole multipolarity. Here 7r M are the conjugate momenta associated to a M . We wish to call 'Bohr 
hamiltonian' the set of more general expressions in which the potential term is a generic function of the parameters. 
We separate the Bohr hamiltonian into three terms: 

H B = T vlb + T ro t + V. (10) 

All the above terms are in general functions of (3 and 7. The first term is the kinetic energy term related with shape 
vibrations with fixed orientation in space, the second instead is the kinetic energy of a rotational motion of the nuclear 
surface without any change of shape. The third term is the restoring potential in the shape parameters. Without 
entering into the details of the derivation [10] we merely restate the definitions of the two kinetic terms using the /3, 7 
variables: 

H 2 (l d a4 d 1 l__d_ . d_\ 

mb ~ ~2B^\J±W c^ + /^sh^a7" 7 97| 



(11) 



3 



Y2 



Trot - ^-^^7 \ T2- (I 2 ) 

fe=1 loir, f„, _ 2tt J 



[sin( 7 -f fc) 



The vibrational part is usually divided into a (3 and a 7 part (corresponding to the two terms in parenthesis in (11)), 
though the two variables are actually mixed by the second term. We will consider potential energies that are only 
function of the internal coordinates (3 and 7. The wave functions ^{(3, 7, 0i) are solutions of the eigenvalue equation 
H B ^ = -E* for the hamiltonian (10). 

It may be thought that, in some sense essentially, the equation to solve is nothing but the Schodinger equation 
and that the record of cases that have already been known in the quantum mechanical context, apply to the present 
situation as well. This is partly true, but there are anyway a number of major differences that must be underlined 
and for which the collective hamiltonian deserves special cares: the Bohr equation is richer being expressed in terms 
of two variables (read two out of five), its 'natural' space is five dimensional instead of three dimensional and these 
affects not only asymptotic behaviours and boundary conditions (that would be trivial), but also the group structure 
that we can identify with the Bohr hamiltonian. 

The main part of the following will be concerned with 7— unstable solution, that is to say, exact solutions with a 
potential that is independent of 7. At the end of this enumeration we will treat approximate solutions with potentials 
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that have also some dependence on 7. 



III. 7 UNSTABLE CASES 

Whenever the potential energy is only a function of (3 the hamiltonian is separable [11]. Setting 

*O9,7,0<) = f(0Ml,Oi) (13) 
we may write two equations, one for the (3 variable, 

that may alternatively be expressed in the canonical form as 

( H 2 / d* ( T+ l)( T + 2) 
\2B m \ 5/3 2 + 02 

and the other for the 7 variable 



(/3 2 /(/3)) = A(/3 2 /(/3)) (15) 



1 5 . _ a q 2 



sm 



sm 37 07 07 4 f— ^ (sm (7 — ^fc)) J 

where A = r(r + 3) is the separation parameter with r = 0,1, 2, ... . Since the potential is only a function of /3, the 

- 2 

angular momentum Q and its third component Q z are constants of motion and the quantum numbers associated 
with them, L and M, are thus good quantum numbers. The so-called 7— angular part of the wavefunction may be 
written as 

*(7,*<) = £ 9 L K Tjy v L MK (e t ) 7 (17) 

K=-L 

where the rotation functions are eigenfunctions of the two operators Q z and Q3, respectively the third component of 
the angular momentum vector along the z axis in the fixed frame of reference (eigenvalue M) and the third component 
of the angular momentum vector along the z' axis in the intrinsic frame of reference of the nucleus (eigenvalue K). 
The functions g, that have the property gK = g~K, are given explicitely by Bes [12] and in addition to L and M, 
are also labeled by other two quantum numbers r and v. The former is the quantum number that comes from the 
solution of the eigenvalue equation of the Casimir operator of the SO(5) group (also called SO(5) seniority), while the 
latter is an empirical label needed in order to distinguish between the occurrence of multiple set of the same L, M 
within a given SO (5) IR and takes the values [13] : 

v = 0,1,2, ... ,[r/3], (18) 

where square brackets indicates the integer part. Here we cannot omit to underline a fact that is worthy of remark: 
Bohr correctly denoted with r the set of these two quantum numbers, but sometimes in the literature the r quantum 
number has been taken as the label of SO(5) (as we do) and the v quantum number has been left out (or forgot!). 
We think that the reason of this omission is that the importance of v is seen only when one takes into account states 
with t > 6. In appendix A we give a detailed discussion of the determination of the sequence of repetitions. 
The K quantum number takes the value r — iv and for each K one has the following list of possible Us: 

L = 2A, 2K - 2, IK - 3, ■ • ■ , K (19) 

that is to say all the integers between AT and 2 A' except for 2A" — 1. The resulting spectrum displays a typical 
degeneracy pattern in energy that is summarized in table I. 

The sets of quantum numbers and the $ eigenfunctions of (16) discussed above are common to all 7— unstable 
problems of this section. The main differences in the spectra are thus to be searched in eq. (15), where in principle 
every possible potential function may be chosen. It is our aim to give here a discussion, as complete as possible, of 
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TABLE I: Quantum numbers for the 7— unstable states. Every group of states with the same r has the same energy. The 
lowest possible multiple occurrence of states with the same L within a given IR of SO(5) is marked in red (See appendix for 
details) . 



FIG. 2: Plot V(/3, 7) of the harmonic oscillator potential (20) discussed by Bohr with C = 1 and < 7 < 3n/2. 



the known cases. 



A. Bohr's harmonic oscillator solution 



The solution given by A.Bohr [1] was historically the first one. He put from the beginning an oscillator potential 
of the form 

V W = \c E I «m I 2 = (20) 

in eq. (15), obtaining the harmonic oscillator hamiltonian in a five dimensional space. The potential V(/3, 7) is 
depicted in fig. 2. The spectrum is given by 

E=(N + b/2)hw (21) 
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TABLE II: Quantum numbers for the 7— unstable states of the harmonic oscillator potential. Every group of states with the 
same N has the same energy. The first repetitions are marked in red. See appendix A and text for details on the sequence of 
repetitions. 



with to = \J C / B and 

N = n A» = 2n /3 + T (22) 

with TV = 0,1,2,3,.... The numbers represent the number of phonons with a given fih component of angular 
momentum along the z axis. The last relation evidences the fact that the 7— unstable case, with a harmonic oscillator 
spectrum and a minimum in j3 = 0, has a further degeneracy with respect to the case discussed above: to a given ./V 
may correspond different sets of (np, r). 



With the content of Table I in mind, we list in Table II the quantum numbers of the first few states. States with the 
same ./V have the same energy. We perform also in this case the analysis of the number of repetitions, that can be 
found in appendix A. 

It is customary to normalize the energies in such a way that the lowest state is at 0, and the first excited is at 1. 
This normalization is equivalent to set the overall energy scale and the relative energy scale respectively. It turns out, 
obviously, that the energy of the second group of excited states in this energy scale is at 2. The ratio R4/2 = -E4+ /E 2 + , 
or better (E A + — E + ) / (E 2 + — E + ) , is the standard reference point for all the solutions of the collective hamiltonian 
and for the comparison with experimental data. Thus, for the harmonic oscillator, we have E 4 + / E 2 + = 2. The 
indexes are referred to the labeling of the states from the bottom to the top. This is unambiguous for the present case, 
but in the following we will encounter many situations in which the degeneration typical of the harmonic oscillator 
will be removed and we will have two indexes, in order to distinguish between different bands (or families) of states 
and between different states within a given family. 

It is useful to introduce the reduced energy e = E^j^- and reduced potential v{f3) = V(/3)^|?p- in eq. (15). In the 
case we are discussing v{j3) = k(3 2 with k = CB/h 2 and hence 



LJLr^JL- (r + l)( r + 2) 

(3 i df3 P df3 p 



P 4 ^~£+ 1 L X ' + /(/?) = (23) 



[1] Alternatively one can take the standard form 



92 + «- (T + 1 " T + a) -i0»>x(fl = o 



a/? 2 p 2 

with x(/3) = /3 2 /(/3). It is indeed always possible to write a linear second order differential equation in its canonical form. 
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FIG. 3: Spectrum of the harmonic oscillator potential. The energy scale (e') is chosen by fixing fiuj — 1. See text and table II 
for explanations. 



with normalized solution 

-1 1/2 

/n T (/3) = 



r(n + r + 5/2) 



P T Ll + *' 2 {p 2 )e-P 2 / 2 (24) 



that contains associated Laguerre polynomials. 

We display in fig. 3 the lowest energy states for the harmonic oscillator. 

B. Wilets and Jean's solutions 

After a very clear and concise introduction to the subject, Wilets and Jean [11] gave the solution in a couple of 
cases: the infinite square well and the displaced harmonic oscillator. 

1. 'Anharmonic oscillator' 

Their aim was to discuss the addition of anharmonicities to the (3 potential and they took as a limiting case what 
they called, with a somewhat misleading designation, 'anharmonic oscillator'. In reality they were the first to treat 
the case of the infinite square well in the form 

V(0) = \ COUSt ' P < @ w (25) 
[ oo , (3 > (3 W 

where f3 w takes a non-null positive real value. They gave the eigenenergies in terms of zeros of the Bcsscl functions 
and they found that the ratio -R4/2 is 2.20, but they said (cited work): 

it represents the extreme case, not realizable in nature. 

It seems thus that they underestimated the importance of this solution. In fact, only recently, Iachello [6] realized 
that this potential may furnish a good description for the shape phase transition between spherical and 7— unstable 
nuclei, although the real form of the potential at the critical point is, to the leading order, a quartic oscillator: the 
infinite square well approximates very well the behaviour of the f3 potential, being very flat around the origin, and 
displays a qualitative agreement with the quickly rising asymptotic behaviour of the quartic potential at infinity. We 
will discuss later his solution, that provides not only the spectrum, but also eigenfunctions and transition rates. We 
will thus describe the complete solution for this potential as it is given in [6], in the appropriate section. 
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FIG. 4: Plot V(f3, 7) of the displaced harmonic oscillator potential (26) discussed by Wilets and Jean with C = 2 and /3o = 0.6 
(these values are chosen merely for purpose of illustration) and < 7 < 37r/2. 

2. Displaced harmonic oscillator 

The same authors treated the modification of the harmonic potential, called displaced harmonic potential, whose 
expression is given by 

V(J3) = \c((3 -(3 ) 2 . (26) 

It was introduced to describe a situation where the minimum along the (3 direction is located at some fixed value /?o 
(see fig. 4). With the substitutions: 

[Buj E . , 2 , x 

one can recast equation (15) in the following form 

if d 2 ( r + l)( T + 2) , , 2 \ . . , . ,„_, 

Wilets and Jean considered the potential formed by the displaced harmonic oscillator plus the 1 jx 1 term and expanded 
it in Taylor series around the minimum x' obtaining 

v(x) = \ ^ (T+1 ^ T + 2) +(x- xo) 2 ^ = v(x') + ^-{x - x') 2 + 0((x - x') 3 ) ■ (28) 

Neglecting terms of order three the effective potential is approximated by a harmonic oscillator and therefore the 
spectrum may be written as 

e = ( n/3 + l/2)o/ + v(x') , (29) 

with up = 0, 1, 2, .... The spectrum looks like the one in fig. 3 with some proper scaling and shift of the energies. 
The eigenfunctions read 

(p(x) = h n0 ((x - x')V^)e-^- x '> 2 "/ 2 , (30) 

where h denotes Hermite polynomials. This is not an exact solution, but it is worth saying that it is a very good 
approximation for large values of xq. 
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FIG. 5: Plot V(/3, 7) of the Davidson potential (31) discussed analytically by Elliott et al. [15, 16] and algebraically by Rowe 
and Bahri [17] with all the constants taken equal to 1 and < 7 < 3tv/2. 



C. Elliott-Evans-Park's solution for the Davidson potential. 



Davidson [14] introduced a potential of the form ar 2 + b/r 2 for the interaction between the constituents of a diatomic 
molecule. Elliott and collaborators [15, 16] employed (though they were not historically the firsts, see section on the 
Warsaw solution!) a similar form in the context of quadrupole deformations and Rowe and Bahri [17] gave a detailed 
algebraic discussion of this potential both for molecular and nuclear spectra (see section VI). We will confine here to 
the direct analytic solution of the Bohr equation with the potential: 



(31) 



^0 P 

that presents a minimum in (3q and A is related to its steepness. Then the differential equation in the 'radial' variable 



2B n 



_1 JLn 4 — - +A (^-^ 

P A dp' dp p 2 J \p Q p 



>f(P) = Ef(P). 



(32) 



By regrouping the terms coming from the potential V(P) it is possible to rewrite the eigenvalue equation as the 
equation for the harmonic oscillator as follows, 



2 -B»- 



1 8 4 d p(p + 3) 

P^dp 1 ' dp p 2 



+ ^-2^1/08) = £7/03) 



(33) 



where p(p + 3) = t(t + 3) + 2AB m p 2 /ti 2 . The solut ion for the five dimensional oscillator is then known to be given 
in terms of associated Laguerre polynomials 



f(P) 



]] T{n+p + 5/2) 6-P-5/2' 

where we have set b = (H 2 Pq /2AB m ) 1 / 4 . The expression of the spectrum is 

E = hu(2n + p + 5/2) -2A, 



(34) 



(35) 



where the non-integer quantum number p is introduced for convenience. Here n should be interpreted as np and 
hui = b~ 2 . By expressing p as a function of r and by introducing a parameter G we can express the spectrum as (see 
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= 2 -r = 

n = 1 



n = n = n = n = 

G = G = 20 G=100 G^oo 

FIG. 6: Spectrum of the Davidson potential. The energy scale is chosen by fixing hu = 1. Notice that this figure is inspired to 
the corresponding fig. 1 in [17], although we have kept the energy normalization consistent for different values of the parameter 
G and we have added the rightmost case to illustrate the SO(6) limit. 



[17]): 



E = huj(2n + 1 + v/(t + 3/2) 2 + G) -2A, (36) 

where we can forget the energy shift and we can set fku to 1. The spectrum exhibits a number of interesting features: 
for G = we obtain again the harmonic oscillator spectrum with the degeneracy of the (n = 0, r = 2) state with the 
(77 = l 7 r — 0) state and so on. For G — ► 00 the ground state band tends to obey the rule r(r + 3) typical of the 
Wilcts-Jean model. The intermediate situation may be used to describe 7— unstable situations where the + of the 
second band typically lies at higher energy with respect to the 2+ of the ground state band. 

These features are illustrated in fig. 6, where the spectrum of the Davidson potential is shown for a few values of 
the parameter G. 

D. Iachello's infinite square well solution or E(5) 

As previously remarked, Iachcllo [6] brought renewed attention on the square well solution, proposing the square 
well as a convenient substitute for the description of the critical point in the U(5)-SO(6) shape phase transition 
between the vibrator and the 7— unstable rotor. With the potential (25) and setting <£>(/?) = /3 3 / 2 /(/3), z = (3k, 
k = y/e and v = r + 3/2, the (3— part of the Bohr equation may be written as a Bcsscl equation whose solutions 
(regular in the origin) are Bessel J function 



ifi + — 

z 



V , N (t + 3/2) 2 



1 



Lp = (37) 
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with Dirichlct boundary condition (ip((3 w ) = 0). This fact determines the spectrum as a function of the £th zero, 
X£. T , of the J T +3/2( z ) function, namely: 




*C.r R- • (38) 



We display in fig. (8) the lowest part of the spectrum. The wave functions are 

/ £iT (/3) = c ^p-^J T+m (^l) . (39) 

The normalization constant may be found analytically from the normalization condition J d(3 fi^ f 2 {(3) = 1. The 
result, containing an hypergeometric function, reads 

„3+2rn2 



J_ kfy T 0l 1 F 2 (2 + T;7/2 + r,4 + 2r;-fc 2 ) 

,2 



c? 2 4 + 2T r(5/2 + r)r(7/2 + r) (40) 



and, although rather complicated, may furnish an alternative to direct numerical computation. 

This extraordinary simple, but nonetheless very successful, solution has been labeled E(5), as the cuclidian group in 
the five dimensional space of the quadrupole variables. This group label may be interpreted in a very straightforward 
way by noticing that the (3— part eq. (14), and hence eq. (37), may be written as 

[tt 2 + U (/3) -e]/(/3) = (41) 

where are the conjugate momenta with respect to the five quadrupole deformation variables, a^. The hamiltonian 
(41) is invariant with respect to rotations in the five dimensional Hilbcrt space associated with the variables defined 
above, and in the special case of the square well cither u((3) — or u{f3) = oo, thus the hamiltonian in the relevant 
region is written as a function of 7r only (plus a constant). Therefore the hamiltonian is also invariant with respect 
to translations in the five dimensional space. The present solution is in summary invariant with respect to the 
transformations induced by the group E(5) that is the semidirect sum of the five dimensional translation and rotation 
groups. 
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FIG. 8: Spectrum of the five-dimensional infinite square well. Energies are scaled shifting the ground state to zero and fixing 
the lowest excited 2 + state to lie at 1. Transition rates are also reported fixing B(_E2,2i,i —> 0i,o) = 100. Adapted from [6] 




FIG. 9: Plot V(P,y) of the finite well potential (discussed by Caprio [18]), with < 7 < 3vr/2. 



E. Caprio's finite square well solution 

The previous analytical solution is a successful approximation of the physics of an entire class of the so-called 
7— unstable nuclei, based on the infinite well potential. Caprio [18] has investigated whether a more realistic behaviour 
of the potential well, taken as a finite square well, alters this view. The finite square well potential is depicted in fig. 
9 and has the following expression 



V(p) = 



0>0u 



(42) 



where j3 w is the position of the wall or step. Eigcnsolutions may be found: inside the well they are again Bcsscl 
functions, while outside the well, taking into account the correct asymptotic behaviour at infinity, the solution may 
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FIG. 10: Evolution of excitation energies as a function of xq normalized to the first excited state (absolute energies in the 
inset). The threshold between bound and unbound states is indicated with a dashed line. Adapted from [18], courtesy of the 
author. 



be written in terms of modified spherical Bessel functions. The complete solution is 

where vq = (2B m /h 2 )Vo is the reduced potential and e^ iT are the reduced eigenvalues, that are obtained by requiring 
the continuity of the wave function and of its first derivative at the position of the step (3 W . One can define a 
dimcnsionlcss energy variable 

/ ex 1 / 2 

77(e) =1-- (44) 



and the parameter 



(45) 



and substitute into the matching condition to obtain a transcendental equation that must be solved numerically to 
determine the eigenvalues. Unlike the previous case the spectrum, displayed in fig. 10 has a finite number of discrete 
eigenvalues, but only those that are lying close to the threshold are appreciably modified with respect to the infinite 
well solution. For the most part the spectrum does not differ substantially from the E(5) case and the same is true 
for electromagnetic transition rates. 



F. Solutions for the Coulomb-like and Kratzer-like potentials 



The analytic solution of the problem with a Coulomb- like potential (fig. 11, a) and with a Kratzer-like potential 
(fig. 11, b) has been discussed in ref. [19]. For these two potentials it is possible to identify the Bohr's equation with 
the Whittakcr's equation, whose solution is known in terms of Whittakcr's functions. Namely, after having introduced 
reduced energies and potential as in the harmonic oscillator case, we can rewrite eq. (15) in its standard form with 
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FIG. 11: Plot V(/3, 7) of the Coulomb-like (a) and Kratzer-like (b) potentials ([19]) discussed by Fortunato and Vitturi, with 
< 7 < 3tt/2. 



the transformation x(/3) = P f((3). The canonical form for the Bohr equation is thus 



X"(P) + u((3) (T+ /3 3 / 2)2 + ^ }*(/?) = (46) 

and will be used in the following to derive analytic solutions in the two cases cited above. 
The Coulomb-like potential reads: 

uc(0)=~, (47) 

with A > 0. With the substitutions e = —e, x = 2^/e(3, k = A/(2y/e) and /1 = r + 3/2, equation (46) takes the 
Whittakcr's standard form [20]: 

S(*) + {-\ + ± + y^}x{*)=0 (48) 
and its regular solution for negative energies is (as in [20]) expressed in terms of the Whittaker's function Mk^(x) : 

(2/1 + 1) /O / X 

XkA x ) = Afr,iX^e- x ' 2 tFx (fi + 1/2 - k, 2/i + 1; x) . (49) 

When the first parameter /t + 1/2 — k = r + 2 — A/(2y/e) is a negative integer, that we call — £, the Whittaker's 
function reduces to an associated Laguerre polynomial and the wave function may be written as 

x {xU = x ^ e ' x/2 j^^ L{ t } ( x ) ( 5 °) 

where the denominator is a Pochhammer symbol. The condition given above fixes analytically the spectrum 

A 2 /4 , , 

£ ^ ^ WTTTW (51) 

A portion of this spectrum is shown in fig. 12. Its most distinctive features are the position of the (4 + , 2 + ) doublet 
at an energy of 1.35 and the presence of a threshold at 1.8 that corresponds to an infinite quantum number. Note 
also that each group of (£, r) states is degenerate with any (£ + k, t — k) group of states with k = 1, r. 

The Kratzer-like potential may be thought as a modification of the former that may be shaped, adjusting the 
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e' r C = 

FIG. 12: Spectrum of the Coulomb-like potential. The energy scale (e') is chosen by fixing the energy of the first two states 
respectively to and 1. The transition rate for 2q i — * Oq~ has been fixed to 100. Some selected quadrupole transitions are 
shown in the figure for simplicity. 



parameters, in such a way to display a pocket at some fixed point. It has the following expression 



-^--Wf-if I- 



where we have shown two possible ways to parameterize this potential. The first is the easiest to use as it will be clear 
from the formula for the spectrum, while the second is related to the geometrical shape of the potential: T> is the 
depth of pocket and (5q is the position of the minimum. These two sets of parameters are connected by simple relations 
that may be deduced from the equation above. In the following we shall make use of the former notation. Equation 
(46) with the Kratzer potential and the substitutions e = — e, x = 2y/e/3, k = A/(2y/e) and [i 2 = (r + 3/2) 2 + B, 
takes again the Whittaker's standard form. The solutions are the same as abo ve with the new substitutions and the 
same arguments apply for the properties of convergence. Now fj, + 1/2 — k = \J (t + 3/2) 2 + B + 1/2 — A/ (2y/e) must 
be a negative integer, — £. The spectrum is in this case: 

A 2 1 A 

e T ^ = = (53) 

(y/(r + 3/2) 2 + B + 1/2 + 2 

with £ = 0, 1, 2, .. to label different families. In reference [19] the evolution of the spectra with /3o and V is studied 
in detail in order to establish a connection between shape of the potential and spectrum. Note, however, that if we 
use the two parameters A and B, we can realize that A does not play any role in determining the scaled eigenvalues 
and therefore the various spectra depend only on B, ranging from the Coulomb-like limit (B = 0) to the 0(6) limit 
(B — > oo). This is illustrated in fig. (13). The threshold varies with B from 1.8 to infinity and the position of the 
4 + varies from 1.35 to 10/3. This makes the Kratzer-like solution very flexible. The parameter B not only fixes the 
position of the 4q , but also all the other excited states including all the f3— bands. 
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FIG. 13: Evolution of spectra with fixed A — 20 and increasing _B for the two limiting situation (B = 0, Coulomb-like and 
B — > oo, 0(6) limit). The first two bands (£ = 0, 1) are displayed with their lowest states (r = 0, 1, 2, ..). The various substates 
are not displayed for sake of simplicity. Some transition rates are indicated. 




FIG. 14: Plot V(/3, 7) with u + (j3) of the sextic potential discussed by Levai and Arias [21] with < 7 < 2>n/2 for r = 0, M = 1 
and a set of parameters chosen for the purpose of illustration. 



G. Levai and Arias' sextic potential solution 



J. M. Arias and G. Levai [21] proposed the sextic oscillator as an example of quasi-exactly solvable potential for the 
Bohr hamiltonian. Only a certain number of eigenvalues may be obtained in closed form for this class of potentials. 
In the present case this is possible for the lowest few values of the principal quantum number. The sextic oscillator 
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TABLE III: Summary of the lowest eigenvalues of the sextic oscillator. For fixed parameters a, b the energies for which a close 
form may be obtained are labeled by n, the number of nodes of the radial wavefunction, and r. 



has the following expression 



u^ip) = (6 2 - 4ac 7r )/3 2 + 2a6/3 4 + a 2 /3 6 + w£ , (54) 

where the parameters a, 6 are used to determine, with an ample choice, the shape of the potential. The index ir = ± is 
used to distinguish between potentials for even and odd values of r. The constants c K are obtained setting r = 2s — 5/2 
(that must be a non-negative integer) and requiring that 

s + M + — = c = const. (55) 



being M a non-negative integer. For each M it is possible to get analytically M + l of its solutions. The last condition 
is needed in order to keep constant the quadratic term of the potential, otherwise the shape of the potential will change 
with s. This implies also that when M is increased by one unity, the value of r must correspondingly decrease of 
two units. The constants Uq may be used to shift the minima of the two potentials at the same energy. There are 
various considerations that can be made here, for which the reader is referred to the cited work. For our purpose it 
will suffice to summarize the solutions in table III, where the formulae for the lowest eigenvalues are given. 



We give in Fig. (14) an example of the shape of the potential surface in order to help the reader to visualize the sextic 
potential. This potential is a rather flexible one, since the parameters may be adjusted in order to yield a minimum 
at j3 = or at /3 > 0. In the latter case it might also have a local maximum at some (3 > before reaching the global 
minimum. These features may be very useful for an accurate description of shape phase transitions. 

The unnormalized solutions of the canonical form of the differential equation (having set <f>(/3) = (3 2 /(/?)) may be 
written as 

MP) - P n {(3 2 )f»-V 2 exp(-^- - M!) , (56) 
where P n are polynomials of order n. The property of normalizability imposes a > 0. 

H. Linear potential 

We give in this paragraph an approximated treatment d la Wilets et Jean that may be applied to the linear potential 

u(j3) = A(3 (57) 

being A a positive constant. This is depicted in fig. 15 

As far as we know this potential has not been treated in any paper concerning the subject. Regrouping all the 
terms except the energy in the second term of eq. (46), we can define an effective potential 

y^m = ^+ {z± ^'^ (58) 

that can be expanded as a Taylor series around the minimum that is located at /3 = (2(r + l)(r + 2)/A) 1,/3 . 
Neglecting terms of cubic or higher orders, we are now in the conditions to solve the differential equation (analogous 
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FIG. 15: Plot V(/3,7) of the linear potential with < 7 < 3tt/2. 
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TABLE IV: Comparison between approximated eigenvalues calculated with the method explained in the text, and the numerical 
values courtesy of M.Caprio [22]. The result is satisfactory for the first family (1 to 3 % errors), while it gets worse, as expected, 
when n growths. 

to the displaced harmonic oscillator) 

x "( ( 3) + {e-^/3 -M(/3-/3 ) 2 }x(/3) = (59) 

that gives the following spectrum: 

3 1 / 3 A 

e =-Apo + (nf3 + -)J—. (60) 
2 2 I ft 

The wave functions are expressed in terms of Hermite polynomials H of order np 

X {fi) = H nf ,((f3-p Q )i/3A^)e-^-M 2 / 2 V^ (61) 

When one measures, as usual, the energies starting from the ground state in units of the energy of the first 2+, the 
spectrum does not depend on the A parameter. 

In table IV we report a list of the lowest eigenvalues calculated in the present approach (by the author) and 
numerically by Caprio [22]. He had performed numerical calculations to check the validity of this approximation as a 
premise to the so-called sloped wall potential, that is flat from the center to a given point after which it has a linear 
dependence (see Sect. (Ill 1 5)). 
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I. Other cases 



The cases treated in the previous pages do not exhaust completely the list of known (analytic or approximated) 
solution, but, to remain loyal to the statements made in the introduction, we preferred to distinguish between the 
former cases and the present subsection. Solutions of the Bohr hamiltonian have been derived, for instance, from the 
classical limit of the interacting boson model (IBM). Other potentials, that are interesting for a phenomenological 
modelization, have been treated numerically. We will briefly examine here some of these cases. 



1. Ginocchio's anharmonic potential solution 

The idea in Ginocchio's paper [24], that we will resume without going into the details because it would be a task 
beyond our aims, is to bridge the IBM and the collective model by means of coherent boson states. Starting from an 
IBM hamiltonian with anharmonicitics (and with a fixed number of bosons) and pairwise interactions he derived a 
second order differential operator which has a spectrum that contains as the lowest levels the IBM cigenenergics. The 
radial profile of the potential connected with this solution is very interesting since it is negative in the origin and is 
very flat near it; it grows rather steadily towards zero at some given point and has a small tail that asymptotically 
approaches zero. It may be written as 

W) = ^IW' (62) 

where the strength of the potential is connected to the number of bosons and b measures the anharmonicity of the 
system. The eigenfunctions are given in terms of Jacobi polynomials: 

<P*M = N nT (l - x) T/2 (l + x) 1+a '' 2 P% a \x) , (63) 
where in = (n — t)/2, a = t + 3/2, a' = N — n and 

*=^£. (64) 



1 + 6/3 2 



The normalization constant is given by 



2 _ Wm!r(a'|tt + m|l) 

n,r ~ 2 «+«'-i/2r(a + TO +l)r(a' + m+l) ' { ' 

The reader is referred to the original paper and to Rcf. [25] for details on the normalization, on the volume clement, 
on the eigenenergies, on the additional definitions and on the quantum numbers that find their roots in the IBM. 



2. Warsaw solution 



In rcf. [26] an interesting 7— unstable case is solved numerically and compared with data. In particular the so-called 
Myers- Swiatccki potential: 

V(f3 2 ) = ^V c f3 2 + G[er^l a ? - 1] (66) 

is solved. With this potential the authors say that they get practically the same results that can be obtained with 
a modification of the Davidson potential. They also give the analytic solution in this case (about ten years before 
the work by Elliott and collaborators). This correspondence can be understood, developing the exponential function, 
from the interplay between the f3 2 and /3~ 2 terms. 

In ref. [27] the generalized Bohr hamiltonian is employed (we do not go into details, see [28]) with a potential of 
the following form: 



V(P) = \c 2 p 2 + C 8 /? 8 + G[e-^/ Q ) 2 - 1] 



(67) 
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4.92 
6.26 



8.01 
9.50 



11.34 
12.95 



r 



3 



TABLE V: Excitation energies for a /3 4 potential relative to the energy of the first excited state (from [29]). 



It is interesting to note that the comparison with the spectrum of 134 Ba that they obtained is not less valuable that 
the one obtained in the context of the critical point symmetry E(5). We would like to point out that, if (again) a 
series expansion for the transcendental function is used, the result is that we have to deal with a generalization of the 
Davidson potential with even powers of (3. This has been treated in [31] where the authors show how the solution 
of such a potential tend to the infinite square well solution when the leading power increases. They have reasons to 
believe that a term like /3 s is enough to have a good correspondence. 



The quartic oscillator in (3 is expected to play an important role in the shape phase transition from harmonic 
oscillator, U(5), to 7— unstable cases , 0(6). In fact it can be shown [29] that the large N limit of the IBM at the 
critical point and the (numerical) solution of the Bohr differential equation with a /3 potential lead to the same 
results and these results are fundamentally different from the analytic solution obtained from the solution of the Bohr 
equation with an infinite square well, the so called E(5) symmetry. Thus it should be stressed that E(5) symmetry 
is the exact mathematical solution for an approximate physical problem. In table V we give numerical values for 
the excitation energies of the /3 4 potential and we refer to [29] for a thorough comparison with E(5). They obtain a 
geometrical limit of IBM by using a coherent state formalism. Considering two different IBM hamiltonian they obtain 
energy surfaces showing that a (3 A potential is associated with the critical point. It is also argued that the IBM could 
provide the finite N correction to the predictions of the simple collective model. This is of interest for an identification 
of nuclei at the critical point. 



A study that deserves comments is a generalization of the harmonic oscillator and quartic potentials that has been 
discussed in [30]. The authors solve numerically the Bohr hamiltonian with the sequence of potential f3 2n with n 
positive non-null integer. 

The harmonic oscillator case occurs for n = 1, while the square well may be considered as the limit for n — > 00. 
The structure of the spectra of all these potentials changes smoothly from one extreme to the other, thus bridging 
the two solutions. This model allows thus for a large variety of different spectra with the R4/2 ratio varying from 2 



The E(5) and X(5) solutions generated not only an experimental effort aimed at the recognition of new patterns 
in nuclear spectra, but also a theoretical effort aimed at the assessment of its most important characteristics. The 
solution for the finite square well and the solution for the sloped wall potential are parts of this effort. In the former 
case the conclusion was that the results of the E(5) model are quite robust with respect to the introduction of a finite 
depth. In the latter case a flat potential with a sloped wall was investigated [32] to establish if the results of the E(5) 
and X(5) symmetries are robust with respect to (small) inclinations of the wall of the infinite potential well. The 
sloped-wall potential reads: 



3. Quartic potential 



4- Bonatsos' et al. solution 



to 2.20. 



5. Caprio 's sloped wall potential 




(68) 



23 



The solution in the internal region may be found analytically in terms of Besscl function, while in the external region 
an analytic treatment is not possible. The author notice that when the term with a /3~ 2 dependence is not present, the 
equation in f3 reduces to the Airy equation for which the solution is possible. The problem is then solved numerically 
using the Airy functions as an efficient basis for diagonalization. The eigenvalues are determined by the condition of 
continuity of the logarithmic derivative at the matching point. The eigenvalues are lowered relative to the respective 
infinite square well cases, and the lowering is more effective for high lying states. On the whole the spectrum is 
very sensitive to the stiffness of the wall, and the deviation from the X(5) predictions are able to generate a closer 
agreement with nuclear spectra in the N ~ 90 region. 

IV. AXIAL 7 STABLE CASES. 

The present section deals mainly with potentials of the type w(/3, 7) = u{(3) + 17(7), where the term in 7 is taken as 
an harmonic oscillator. This clearly violates the property of periodicity in 7 that our problem possesses. However if we 
consider only a narrow interval of 7's around zero it is possible to give an analytic solution to the differential equations 
using the expansion of the periodic functions around 7 = (see for example [11] for a thorough analysis). Although 
this is not the only possibility to obtain solvable, or approximativcly solvable, Bohr hamiltonians, we will discuss it 
in detail for its importance in connection with the issue of critical point symmetries and shape phase transitions. 
This was first discussed in ref. [7] , where the approximations that will be used throughout were introduced. Within 
this class of potentials one may study approximate solutions for a /3— soft, 7— soft rotor. In subsection E a different 
method, that leads to an exact separation, is used. 

It is perhaps worth saying that when the 7 degree of freedom comes into play, one usually have to consider 7 
phonons. The f3 phonon shares the same total angular momentum (A = 2) with the 7 phonon, but while the former 
has a null component along the quantization axis, the latter has a non null (^1 = ±2) component. The rules to assign 
quantum numbers to the various states, for a given number of quanta in 7, are: 

j = 0, ...,n 7 
K = 2n 7 - 4j 

r 0,2,4,6,8,... K =o 

\ K,K + l,K + 2,... K ^0 
When K is negative one must consider its absolute value. 

A. Iachello's square well solution or X(5) 

This important case gave the starting signal to a number of experimental works about the so-called X(5) symmetry. 
Iachello proposed an approximate separation of variables for the Bohr hamiltonian that consists of two steps. The 
first point is to realize that around 7 = 0° the rotational kinetic energy, eq. (12), may be written as 

Y ^ 2-i(Qi + Q 2 2 + Ql) + Ql(^-t)- ( 7 °) 

^sin( 7 - ffc) 2 ^Z1_?X_^ > ^ 3 Vsm 7 2 i) 

In this case the problem is no longer 7— unstable and the separation of the wavefunction d la Bes do not hold anymore. 
We should instead look for solutions of the type ^>((3,j,8i) = <Pk(P, j)T>m _r-(^)j where V is a Wigner function, 

eigenfunction of the square of total angular momentum and of its third component. The action of L 2 and Q\ leaves 
the following equation 



[2] This subsection has been inserted here because of its numerical character, although a better place would have been at the end of the 
next section. 
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J__d_ ,d_ 1 l_d_ . d 1 



3 V sin 7^ 3 



+ U (/3, 7 )^|(/3,7) = 7). (71) 



Considering now a potential of the type 7) = it(/3) + 11(7) the above equation may be approximately separated 
into the following set 



1 d 4 d L(L+1) 

/3 4 d(3 P 8/3 3/3 2 



■«(/?) 



£l(/3) = e/3a(/3) 



(72) 



(/3 2 ) sin 37 c?7 



9 . „ d 
sm 37— 



9 7 4(/3 2 ) 



A'" 



sin 7^ 



1^7) 



(73) 



where (/3 2 ) is the average of (3 2 over £(/?). Here we have e ~ e/3 + e 7 and yff(/3, 7) — £l(0)tik (7)- Some authors 
noticed that a different set of equations, in which the term K 2 / (3(f3 2 )) is kept in the first equation, may have been 
derived as well. In this case, that we will not take into account, the solution of the (3 equation would carry also the 
K quantum number. 

Insofar the procedure has been carried for a general potential. Iachello treated a square well potential in /3 combined 
with a harmonic oscillator in 7. Introducing the square well potential (25) into cq. (72) and setting = /3 3/2 £(/3), 
ep = k 2 p and z = /3kp, he obtained a Bessel equation 



e + - 

z 



'-4 

z 1 



where 



1) 



9\V 2 
4/ 



(74) 



(75) 



The solution of this equation is therefore expressed in terms of Bessel function with a non-integer index, v, 

6,l = c s , L p-V 2 J u (k s , L /3), (76) 

where the set of quantum numbers refers to the sth zero of the Bessel function J v (z) and to the total angular 
momentum, L. In fact the boundary condition at the wall of the well, £(/3 w ) = 0, determines the eigenvalues, 
because the wave function has a node if and only if the Bessel function has a zero, that we call x s> l- The spectrum 
is thus 



e s ,L 



(77) 



where k Sj i = and (3 W is the position of the wall. 

We give here the solution of the 7— part for the harmonic oscillator, as it was given by Iachello. The following 
subsections IV B - IV D will share the same treatment for the 7— variable. A trigonometric expansion, strictly valid 
in a small sector around the origin (sin 7 ~ 7), is made in eq. (73) with a harmonic oscillator potential, obtaining 
formally the radial equation of a two dimensional harmonic oscillator: 



1 1 
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(78) 
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The solution is given by 
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TP 



= (n 7 + 1) 



(79) 
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FIG. 16: Plot V(/3, 7) of the square well potential in (3 combined with a harmonic oscillator in 7 discussed by Iachello, 
-7r/4 < 7 < 7r/4. 



with 



and 

Vn^Ml) = c„, if 7 |K/2| e- 3 ^ 2/2 ik K| (3a7 2 ), (80) 

where n = (?i 7 — | K |)/2 with n 7 = 0, 1, 2, . . . and iif"' are Laguerre polynomials. The parameter a is the strength 
of the harmonic oscillator potential, or in other words its amplitude. 



B. Solutions for Coulomb- and Kratzer-like potentials 

With the substitution = Xl(P)P~ 2 , equation (72) may be simplified to its standard form: 

w xlW + - u(p) - f — ^ Wm = . (si) 

By inserting the Coulomb-like potential u{[3) = — A/ (3 with A > and recasting the problem in the variable 
x = 2(3y/e, with the further substitutions e = — e, k = and /i 2 = (| + L ^ L 3 +1 - ) ) , we obtain the Whittaker's 
equation 



_d 2 

dx 




!^*Xl(,x) + {— a +- + ^-^}xl{x) = 0. (82) 



Its regular solution is the Whittaker function M^^(x) that reads: 

X l(x) = Me- X l 2 x^ X F Y (1 + ii - k, 1 + 2/i, x) (83) 

which is, in general, a multivalued function. The constant Af is determined from the normalization condition. Wc 
adopt usual conventions and thus the function is analytic on the real axis. The hypergeometric series is an infinite series 
and to recover a good asymptotic behaviour we must require that it terminates, i.e. that it becomes a polynomial. 
This happens when the first argument is a negative integer, —v, that must thence be regarded as an additive quantum 
number. This condition fixes unambiguously the f3 part of the spectrum: 
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FIG. 17: Plot V(/3, 7) of the Kratzer-like potential ([33]) combined with a harmonic oscillator in 7 discussed by Fortunato and 
Vitturi, with — 7r/6 < 7 < it/6. 



Fig. 18 shows the spectrum of the axial rotor with coulomb-like potential. 

The case of the Kratzer-like potential is treated in a very similar way. Inserting in eq. (81) the potential 

uifi) = -A/ /3 + B/[3 2 



(85) 



and setting x = 2/3-^/e, e — — e, k = and /i 2 = (B + | + L 3"^ )) we obtain again the Whittaker's equation 
whose solutions may be written as in eq. (83) and we may repeat the whole procedure of the previous section, the 
only major difference being the definition of fj,. The spectrum assumes in this case the form 



E (/3) 



A 2 /A 




(86) 



L(L+1) 1 
3 '2 



We recall (see sect. IIIF) that the two parameters used here, A and B, have an immediate translation into the 
position, /3o, and depth of the minimum of the potential, T>, being valid the following relations: A = 2(3qT> and 
B = ffiV. Consequently /3o = 2B/A and V = A 2 /(4B). 

In fig. 19 we report a study of the evolution of the spectrum and transition rates for the Kratzer-like rotor as a 
function of the parameter B. When B is large a typical rotational spectrum is recovered. 
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FIG. 18: /3— spectrum of the Coulomb-like case (n 7 = 0). The vertical scale is expressed in units of the energy of the 
(v — 0, L = 2) state. Some B(E2) values are displayed in units of the lowest transition of the first band (numbers in italics). 
Notice the presence and energy of a threshold and the absence of degeneracy. From [33] . 



C. Bonatsos' et al. solution 



As in sec. Ill 1 4, we summarize here briefly another work by Bonatsos and collaborators [31] that deals with the 
numerical solution of a sequence of potentials interpolating between the U(5) and X(5) models of the form: 

«2»G9,7) = —+ c ^> ( 87 ) 

where the harmonic dependence on 7 is the same of the X(5) case, while the potential in may be considered as a 
generalization of other cases. For n = 1 one gets an exactly soluble model that is called X(5)— 2 . while for n — > 00 
the infinite square well potential is recovered. As expected energy ratios and transition rated change smoothly from 
one limiting case to the other and it is shown that the X(5) results are already well approximated for n ~ 4. 



D. Pietralla and Gorbachenko's solution or CBS-model 



An extension of the X(5) model has been recently proposed [34] with the aim to study the evolution of the first 
excited + state in transitional nuclei (N ~ 90). The authors consider the same problem discussed by Iachello for 
axially symmetric prolate nuclei, but they take as a potential in f3 the infinite square well with two boundaries: 

{OO , /3 < 0m 

, m < 13 < I3 M (88) 
OO , 0>0M 
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FIG. 19: Evolution of the f3— spectrum of the Kratzer-like rotor with the parameter B. The vertical scale is expressed in units 
of the energy of the (v = 0, L = 2) state. The B(E2) values, calculated with formula (??), between the lowest states are 
indicated in italics beside the corresponding downward arrows. Intraband transitions are not marked for the sake of simplicity. 
From [33]. 



FIG. 20: Plot V(f3, 7) of the two-sided infinite square well potential in j3 combined with a harmonic oscillator in 7 discussed 
by Pietralla and Gorbachenko [34], with — rr/12 < 7 < 7r/6. 



The number rp = (5 m j (5m may be used to parameterize the stiffness of the potential. The X(5) solution is a special 
case of this model for rp = as well as the rigid rotor is another special case for rp = 1. Eq. (72) may be transformed 
into the Bessel equation with the substitutions z = yj Ej (h 2 /2B m )j3 and £(z) = /3 3 / 2 £i(/3). The solution is a linear 
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combination of Bcsscl's J v (z) and Y„(z) functions with index v = \J L(L + l)/3 + 9/4. The authors give the complete 
solution in the following form: 

= jjjjfc [M^lJIPm) + 7YY u (z r j? a p/p M j\ , (89) 

where cj, jS is a normalization constant. The spectrum is found imposing the condition that the wave functions in p 
must have nodes at the two boundaries: 

e(/?m) = Z(Pm) = o 

or 

i(rpZm) = i(zM)=0, (90) 



where Zm = \J E/(h 2 /2B m )P^[. The quantization condition is thus 

J v {zm)Y v {t0Zm) - Jv{t0Zm)Y v {zm) = , (91) 
whose sth zero, z^ a , must be calculated numerically. The eigenvalues are found as 

El, s = ^rizZ) 2 . (92) 

This solution is rather interesting because there is one more degree of freedom with respect to X(5) and moving rp 
from to 1 has not only the effect of increasing the position of the P bands, but also to invert the trend of the ratio 
i? 4 / 2 (s): while for rp — the ratio R4/2 of the ground state band is higher than for any excited band, for greater 
values of rp the opposite happens. Another interesting feature of this model is the possibility to study the effect of 
the 'centrifugal stretching': the term with P~ 2 power acts as a 'centrifugal stretch' (although for the sake of clarity 
the variable here is not the radius of the system, but its deformation) in the sense that the wave functions of states 
labeled by high values of the SO (6) quantum number will be pushed towards higher deformations, decreasing their 
probability amplitudes in the region close to the origin. Within this model one actually cut the region close to the 
origin with the internal wall. As a consequence only states with appreciable amplitude in the classically forbidden 
area will be more affected: the overall result is that states with high r arc not much affected by the presence of 
the internal wall, while states with low r arc shifted up in energy. The spectrum will therefore show smaller energy 
distances among low-lying states with respect to X(5). 

The authors named their model as "confined /3— soft" rotor model (CBS) by analogy with the Wilets-Jean 7— soft 
model (see Appendix B). 



E. Prolate and oblate axial rotor: new cases 



We propose in this section new solutions that, as far as we know, have not been published anywhere else, although 
they are pretty simple. To obtain them we have done nothing but collecting knowledge coming from the various 
solutions that have been devised so far and apply known approximations. In doing so we made use of some suggestions 
due to D.J.Rowe [35]. The novelty here is represented by the transformation (97) that allows a better level of 
approximation for the goniomctric functions. This is reflected in four new analytically solvable cases. 

The Schrodinger equation for the Bohr hamiltonian may be profitably simplified for the prolate and oblate soft axial 
rotors (respectively around 7 ~ and around 7 ~ 7r) using the approximation (70) devised by Iachcllo. Alternatively 
one can treat the 7 ~ 7r/3 case by noticing that in that case the projection of the angular momentum on the intrinsic 
x-axis is a good quantum number (see fig. 1). We will seek solutions of the type 

*G9,70i) - tdP^Kd^M^m . (93) 

Since the rotational part is standard and the action of the square of the angular momentum and of its third component 
are trivially evaluated, we will restrict to the P — 7 part, as done in [7]. The result is formula (71), that we employ 
here as a starting point for new solutions. 

Contrary to what was done in the X(5) model [7], we may separate exactly the Bohr equation as explained in 
section VD, using a potential of the form u(/3, 7) = u\(p) + it2(7)//3 2 , obtaining the following two equations 
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1 d d K 2 

■ sin ( 3 7)— + — —2 ~u >Vk(i) = (95) 



sin (37) ^7 97 4 sin (7) 

where u is a separation constant and the energy, e, is contained only in the first equation. Notice that one could have 
alternatively retained one or both of the terms L(L + l)/3 and K 2 /3 in the equation for 7. We prefer to keep them 
in the (3— part because this part is solved exactly. The 7— part usually experiences a second turn of approximations. 

In the second equation one can adopt the same trigonometric simplifications (sin 7 ~ 7 and cos 7 ~ 1) that are 
implicitly used in (70) or may try a more sophisticated approach. In fact eq. (95) may be written as 



d 2 n n d . . K 

_+3cot 37^+0,-^(7)-^ j (7) 



2 +3cot3 7 — +lo-u 2 { 1 )- W7)=0 (96) 



and with the transformation 



Vk(V = A1 — q- (97) 
ysin 37 

it may be brought in the standard form 

+ "> ~ « W + I ~ 7^ + i cot2 (37)}pk(7) = , (98) 

where the term with the first derivative has been eliminated. It is worth noticing that now a better approximation 
may be used for the trigonometric functions, without loosing the possibility to solve exactly the equation. The series 
expansion for the cotangent, including the second order, may be used 

-cot 2 (3 7 )~^-- + — 7 2 . (99) 
4 V " 4 7 2 2 20 ' y J 

Therefore we can recast the differential equation above as 

id 2 27 1 - K 2 1 

} + - U2 ^ + ~^f- \PK{1) = • (100) 



This expression represents a better level of approximation with respect to any equation used so far for the 7 part of 
the problem. The 'extra' term with dependence 7 -2 has the same behaviour of a 'centrifugal' term and the differential 
equation can be solved exactly with harmonic or Davidson potentials for 1/2(7)- 

With 02(7) = C7 2 , defining fx = q 2 /(2\) and A 2 = C - § and t(t + 1) = (K 2 - l)/4 we obtain 

PK{l) = AA 7 t+1 e-^ 2 x Fr (*±^Z£, t + 3 /2; A7 2 ) , (101) 

where the confluent hypcrgcomctric function is not divergent only when it is a polynomial: this happens if the first 
argument is equal to — n 7 with n 7 = 0, 1, 2, 3, .... This requirement sets the formula for the spectrum 



uJ n -,,K = \/C- §J(2t + 3 + 4n 7 ) -3. (102) 



20 

Notice the constraint on the value of C that comes from the square root. With 



(103) 
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FIG. 21: Graphical representation of the values of u) n k for the lowest states of a few bands, obtained from formula (104) with 
A = 1 (the shift, —3, has been omitted). The assignation of quantum numbers to each band follows from (69). 



we can write 



WtH.x = Jc-^(±\K |+2 + 4n 7 )-3. (104) 

The structure of the various bands predicted by this formula is depicted in fig. 21. We shifted the energies eliminating 
the term —3 and we used a constant strength (A) equal to unity for purpose of illustration. 

The expression for u> must be now used inside the equation for the (3 part. Choosing for example the harmonic 
oscillator: u\{(3) = k(3 2 one has: 



e = Vk(2n p + t + 5/2) (105) 

with t is determined through the relation 

, L(L + 1)-K* 

t(t + 3) = Lo n -,,i< H jj (106) 

(one can employ the Davidson potential in a very similar way) . Alternatively we may consider the Coulomb potential, 
U\(j3) = —A//3, or the Kratzer potential, u\{(3) = —A/ (3 + B/[3 2 . In the latter more general case we obtained 

A 2 /4 , 

(107) 



(J\/A + B + oj ni .j { + L{L+ ^- K2 + 1/2 + n p ) ' 



where one should insert the values of w previously found in formula (104). The most relevant implication of this 
model is the difference between bands with opposites signs of K. 



V. TRIAXIAL CASES 



We divided this section from the previous one, although they both deal with 7— stable cases, because in triaxial 
nuclei the minimum of the potential along 7 is not located at n^, with n G Z. This implies the absence of axial 
symmetry. 

In triaxial nuclei the eigenvalue of the projection of the third component of the angular momentum is no more a 
good quantum number (this can be seen already at the classical level). In the special case of 7 = 30° the projection 
of the angular momentum on the intrinsic y-axis, called sometimes R. is a conserved quantity, but in the regions 
0° < 7 < 30° and 30° < 7 < 60° is not possible to use K or R to label bands. Notice that in the axial case the lowest 
member of a group of states with the same L has always the lowest possible value of K, while at 7 = 30° the lowest 
state with a given L has the highest possible value of R (this is a trivial consequence of formula (116)). 
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A. Davydov's classical solution or rigid rotor model 



Traditionally the name of Davydov is associated with a series of papers [36-39] in collaboration with other various 
researchers. Davydov and Filippov proposed the existence of triaxial nuclei, the properties of which have been 
investigated in the adiabatic approximation, which assumes rotation of the nucleus without change of the intrinsic 
state. The equilibrium shape is a triaxial ellipsoid, whose hamiltonian, frequently called the rigid rotor hamiltonian, 
may be written as: 

L «=1 Sin (7 - gTT/t) 

where Q\ are the projections of the total angular momentum operator on the axes of the intrinsic system. The wave 
functions are expanded in terms of (axial) rotational wave functions 



*JM - I JK ) A K , (109) 



A' 



with 



where K = 0, 2, 4, ... and 



(in) 




Within this rigid-rotor approach analytic expressions may be derived for a few rotational levels with low total angular 
momentum, although this model completely neglects the shape vibrational degrees of freedom. 

A more complicated version, containing the kinetic energy in the (3 variable along with the above rigid rotor 
hamiltonian and a potential term V(P) may be employed [38]. The corresponding Schrodingcr equation for this 
7— rigid, /?— soft model is separable and leads to lengthy expressions for the energy of collective nuclear states. Wc 
give here a simplified version of the spectrum: 

Ej, K ,n = (np + 1/2)E P + ci(J(J + 1) - K 2 ) + c 2 K 2 , (112) 

where 71/3 = 0,1,2,.. is the number of (3 phonons. 

Davydov proposed also a solution that takes into account and 7 vibrations. When in the Bohr hamiltonian a 
potential of the type 



= ic(/3-/3 ) 2 + ic 7 /3 2 ( 7 -7o) 2 



is employed, we can look for solutions that are product of f(/3) and $(7, 0j) functions, obeying the following equations: 
( 1 d . „ 8 Q 



sin3 7 |- + . ; ®\ o + £(7 - 7o) 2 - A) $(7, 6,) = (114) 



and 



sin 3 7 <9 7 dj 4 ^ sin ( 7 _ | ^ 



1 nra a ^2 , ^ 2 (A + 2) 



2B m ^ ' 2 W-W + -^--E\p m =* 1 (115) 

where C 7 = lo 2 B contains the frequency of the 7 vibration. 

The above equations can be solved exactly for a spherical nucleus, reproducing known results. Davydov proposed 
also a solution for non-axial nuclei when small oscillations around the equilibrium position are considered. The main 
point is to approximate the effective potential, that consists of the displaced harmonic oscillator and the 'centrifugal' 
term, expanding it as a power series as Wilets and Jean did for the 7— unstable case. In the present approach (as 
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FIG. 22: Spectrum of the rigid triaxial rotor at 7 = 30°, according to the formula by Meyer-ter-Vehn. This spectrum is also 
predicted by the Davydov rigid model, although only a few low-lying states may be calculated analytically. Quantum numbers 
and energies are reported on the left and right respectively. Notice that here the first excited states lies at 3. 

in the previous one, that has only (3 vibrations) a rather complicated expression for the energy levels is found that 
accounts for rotations as well as f3 and 7 vibrations. Although we do not enter into the details, it should be said that 
this approach had considerable success and applications to nuclear spectra. 

B. Meyer-ter-vehn formula 

In a work focused on the spectrum of an odd nucleoli coupled with a rotating triaxial core [40], Meyer-ter-vehn 
introduced also a simple analytic formula for the spectrum of the rigid rotor at 7 = 30°. At that angle two of the 
three moment of inertia are equal (albeit the axis of the ellipsoid does not have equal length) and the hamiltonian 
(108) is axially symmetric around the intrinsic y-axis. The energy spectrum is easily written (apart from a factor) in 
the following form 

Ej,r = J(J + l)-\R 2 (116) 
where R is the sharp projection of the angular momentum on the y-axis. The wave functions 



do not depend on the 7 variable, but only on the Euler angles. 
We illustrate in fig. 22 the spectrum predicted by eq. (116). 

C. Iachello's solution or Y(5) 

The solution called Y(5) was introduced in Ref. [41] with the aim to describe the critical point of the axial to 
triaxial shape phase transition. The starting point is again the set of differential equations (72, 73) obtained inserting 
in the Bohr hamiltonian a potential of the type V(/3, 7) = u(/3) + ^(7) that allows an approximate separation of 
variables. In the present case a displaced harmonic oscillator potential in the variable (3 is considered, namely 



«(/?) = y^-/?o) 2 , 



(118) 
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FIG. 23: Plot V(/3, 7) of the potential used in the Y(5) symmetry discussed by Iachello ([41]): an infinite square well in 7 with 
< 7 < 7uj = 7r/6 combined with a displaced harmonic oscillator in j3. 



while an infinite square well potential is considered in the asymmetry variable 7, as follows 




(119) 



This particular choice of the potential, depicted in fig. (23), is an approximation to a more general expression, 
— cos(37) + £cos(37) 2 , that changes smoothly from an axial to a triaxial minimum, passing through a critical point 
when £ = 1/2. When at the critical point the shape of this potential is flat around 7 = and eventually rises steadily 
at some point: this can be roughly approximated by a square well. Although this is a rather crude approximation, 
it has the advantage of generating an exact solution, that may again be used as a benchmark. When 7 << 60° one 
may also take sin(7) ~ 7 and eq. (73) may be written as 



15 d 

7C*7^c>7 \ l J v 7 

with K = 0, ±2, ±4, .... This equation is the Bessel equation with solutions given in terms of Bessel functions as 



(120) 



Vs,K(l) = CsJ<JK/2{k s ,K~f) 



(121) 



with k s ^K = Xs.kHw Here x St K is the sth zero of the Jk/2 Bessel function. The cigenfunction must be zero at the 
wall (~f w ) of the infinite square well and this fixes the spectrum to depend upon the zeros of the Bessel function as 



^,s,K = k SjK - K /2. 
The full solution is obtained solving also the (3 part of the problem. 



(122) 



D. Solution around 7 ~ n/6 

It has been shown in [42] that whenever the potential is chosen as 

V 2 (l) 



f(/3,7) = W) + 



(123) 



the Schrodingcr equation Hb^((3, 7, 0i) = E^{j3, 7, 6{) is separable [11]. The set of second order differential equations 
that comes from this separation contains the separation constant f2 and reads: 



(124) 
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FIG. 24: Plot V(/3, 7) of the potential used the solution of the 7 ~ n/6 soft triaxial rotor discussed by Fortunato ([42]): a 
harmonic oscillator in 7 with < 7 < 7r/3 combined with a Kratzer potential in j3 (both taken as a very shallow for purpose 
of illustration). 



(-asbisW £ sin w £ + 14M ' " + .5, ..-hlfl )*' »» " (125) 

with ^(/3, 7, 0j) = /(/3)$(7, 0$). The set above may describe a /3— soft, 7— soft triaxial rotor with a potential that has 
a minimum located in /3 = /3o and 7 = 70 = 7r/6. The same considerations about the privileged role of the projection 
of the angular momentum along the intrinsic y-axis apply in the present case. Restricting ourselves to a small region 
around 7r/6, we multiply the two equations above by 2B m /ti 2 and we define reduced energy and potentials. The new 
set of equations reads 

{-*m*m +u M- E+ %) m = * (126) 



("- TTr|- sin ( 3 7)|- + ^(7)-^ + 7 E • 2, Q \ /oO *Mi) = °- ( 12? ) 
V sin (37) 97 c?7 4 snr (7 - 27TK/3) / 

The spectrum is determined by the solution of the first differential equation in which w, that it is found from the 
solution of the second differential equation, plays the role of the coefficient of a 'centrifugal' term and, as we will 
show, yields a non-trivial expression for the energy levels. Around 7r/6, setting 7 = ir/Q + x, the rotational part of 
the Bohr hamiltonian becomes 

E . 2( % m = 4(0? + QI + Ql) +Qi ( — ^ - 4) . (128) 

sin (7 - 2fK/3) v w ' \cos z (a;) J 

" " L 2 

Changing variable in eq. (127), introducing the harmonic dependence of the potential, u 2 (x) = Cx 2 , and using the 
simplifications 

sin 37 = cos 3a; ~ 1 cosx^l (129) 
together with relation (128) leads to a simplified equation 

l^^+Cx 2 -u; + Q-^Q 2 1 U(x,e i ) = 0. (130) 
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The wave function <J>(x,#i) may be written, following [40], as 



V i6 Ja + L o) Vn -" L > R{x) \Pm, r ^) + i-^MUm] (i3i) 

where the angular part is written in terms of Wigner functions labeled by the projection of the total angular momentum 
on the intrinsic y-axis, R, that is a good quantum number, while the functions r\ n ,l,r(%) are eigenfunctions of the 
one dimensional Schrodingcr equation for the harmonic oscillator. The index n 7 is the quantum number associated 
with the vibrations in the 7 degree of freedom. The spectrum is therefore written as 

wi,ji,n 7 - \/C(2n 7 + 1) + L(L + 1) - ^R 2 (132) 

where the first term corresponds to the 7— vibration, while the other two terms correspond to the Meyer-ter-vehn 
formula [40] that accounts for the rotational energy of the 7 = 7r/6 rigid triaxial rotor. By imposing /(/?) = /3~ 2 x(/3)> 
we can write the eigenvalue equation in the (3 variable as 

x"(/3) + (e u x {P) ^)x(P) = • (133) 

The potential ui(/3) is taken of a Kratzer-like form -A//3 + B//3 2 [19, 33]. Setting z = 20y/e, e = -e, k = A/(2y/i) 
and /x 2 = 9/4 + B + u, eq. (133) may be rewritten as the Whittaker's differential equation: 

r d 2 1 k 1/4- u 2 } , , , . 

The reduced eigenvalues are 

A 2 /4 

e(n 7 , np, L, R) = - r — , (135) 

(y/9/4 + B+u LiR , ny + 1/2 + np) 

where w comes from eq. (132). 

Setting as usual the energy of the ground state to zero and the unit of the energy scale to be the energy of the lowest 
excited 2 + state of the ground state band leads to a spectrum that does not depend on A. The interplay between 
the two parameters B and C may be exploited to fit experimental data. Notice that B and C do not separately fix 
the position of the respective (3 and 7 bands, but they both take part in a non-trivial way to determine the energy 
levels. It is clear that, while the solution of the 7— angular part of the problem gives a straightforward extension of 
the rigid rotor formula in which a simple harmonic term for the 7 degree of freedom appears, the full spectrum is 
rather a more complicated function that essentially depends on the choice of the potential in (3. 

We display in fig. 25 and 26 (see captions, adapted from [42]) the behaviour of some energy levels with respect to 
the parameters of the potential. The Davydov (rigid) model is recovered when B — > 00 as can be seen from fig. 25. 



E. Bonatsos' at al. solution or Z(5) 

A critical point symmetry for a shape phase transition from prolate (7 = 0°) to oblate (7 = 60°) shapes was 
introduced in [43] and called Z(5). The authors argue that in such a transition the triaxial region is crossed and the 
middle lies at 7 = 30°. Unless the previous solution the present one does not use an exact separation of variables, 
but rather an approximate one with 7) = u((3) + ^(7), where the (3 potential is taken to be a square well and 
the 7 potential is taken as a harmonic oscillator with the minimum in 7 = 30°. The infinite potential well in (3 may 
be thought as the critical point between a triaxial vibrator (with the minimum close to (3 = 0) and a triaxial rotator 
with a minimum for a non null value of (3. In this aspect the model closely resembles the X(5) model. 

The same trick adopted in the previous section (formula (128)) is used here to treat the rotational degrees of 
freedom. The Schrodinger equation is approximatively separated into the set 
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FIG. 25: Reduced energies of the lowest state of the /3— band (dashed line) and of a few lowest states of the ground state band 
(solid lines) as a function of B. The limits for the energy levels when B — » oo, that correspond to the rigid triaxial rotor 
energies, are reported in the right side. Here we fixed C = 1. 
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FIG. 26: Reduced energies of the (J* = 2 4 



l,H/9 = 0) state (solid line) as a function of B, for various values of the 



strength of the harmonic potential in 7, C. The first 2 + of the ground state band (dashed line) is reported for reference. 



and 
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— — -— sin 37— + u(7) 77(7) 
sin 37 07 07 
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where the angular momentum quantum number. L, and its projection on the intrinsic y-axis, R, are explicitely 
contained in the first equation. In the second equation the average of f3 2 over appears: this approximation is 
strictly valid only when the potential u{(3) is deep and consequently the mean square value of (3 does not oscillate 
much and remains constant also for a few lower excited states. The energy is e ~ ep + e 7 . With the transformation 
£(/3) = /3 3 / 2 £(/3) and the substitutions ep = k"j } and z = fikp, eq. (136) becomes a Bessel equation with eigenfunctions: 



(138) 
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FIG. 27: Plot V{p, 7) of the potential used the so-called Z(5) solution discussed by Bonatsos and collaborators ([43]): a 
harmonic oscillator in 7 with a minimum in 70 = n/6 combined with a square well potential in /3 (the lowest corner of the 
frame box corresponds to the origin of polar coordinates). 



where the non-integer index is 



y/AL(L+l)-ZR 2 + % 



(139) 



The boundary condition at the wall of the potential well, £(/3 w ) = 0, determines the spectrum as a function of a; Sj „, 
the sth zero of the Bessel function J v {z) 



= i^-f) ■ (140) 
The 7— part is obtained, limiting ourselves to small oscillations around the minimum, from the solution of the equation 

(141) 

where 7 = 7r/6 + x. The above equation is an harmonic oscillator equation with eigenvalues 



2c . 

ir(«7 + l/2), 



7 v & 

with n-y = 0, 1, 2, ... . The eigenfunctions are expressed in terms of Hermite polynomials 



(142) 
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H n (bx)e- b2 * 2 / 2 



(143) 



with b = (c(/3 2 )/2) 1 ^ 4 . The index of the Hermite polynomials is now the number of 7— phonons. 



F. Jolos' solution 



One of the main objectives of ref. [5] is to summarize the various researches on shape phase transitions in nuclei. 
Using the interacting boson model and the coherent state formalism the author analyzes the phase diagram of a cold 
nucleus discussing the order of phase transitions. Within the same perspective the author summarizes also some of 
the recent achievements in the solution of the Bohr hamiltonian. In particular he proposes an approximate solution 
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for the critical point of the phase transition between spherical and triaxially deformed shapes. A potential of the form 

(144) 



F(/3, 7 ) = U (/3) + i^ 6 cos(3 7 ) 2 



is used. Assuming D large enough, one can consider only small oscillations around 7 = it/ 6, and defining x = 7 — 7r/6, 
the equation becomes 



2B n 



J_ d_ A d_ J_(d^_ 9DB, 
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m /3 s x 2 - L(L + 1) + -R 2 ) + u((3) - E W(/?, 7 , 9$ = . 



(145) 



where R is once again the projection of L on the intrinsic y-axis. If the wave function is factorized in the following 
way 
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(146) 



where uj = 2(D/B m ) 1 / 2 and H n are Hermite polynomials, (having used the original notation), one can exploit the 
action of the rotational and differential operators in the variable x. The equation for the (3— part of the problem is 
therefore 



2B n 



f3 4 d(3 P dp /3 2 



L(L + l) 
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(147) 



where one can use the infinite square well model for the potential u(/3) and take the assumption /3 2 — > ((3 2 ). Hence 
the wavefunction for the (3 part may be written as 

/(/?) = /r 3 /V iiR j„(4y?), (148) 

where (notice the slight change of notation with respect to the original) 



L{L + 1)-\r> + \ 
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with (3 W is the position of the infinite wall and x\ R is the i— th zero, i = 1, 2, 3, ... of J v . Here c\ R = V2/p w J' v (xi R ). 
The energy is 



E 



L.R 



2B m pl 



{xl r) + hujn(iLRn \ P 2 \ iLRri). 



(151) 



VI. ALGEBRAIC METHODS 



We have preferred to dedicate a separate section to the solutions of the Bohr hamiltonian obtained through group 
theoretical techniques, even if this breaks the chronological order and necessarily implies a redoubling of some topics, 
because of the insight that they provide. We focus here especially on the results obtained with the su(l,l) algebra. 
More informations on this argument may be found in refs. [44], ch. 18 and [45, 46], while a thorough group theoretical 
analysis of the collective model is found in ref. [47, 48]: the works of Chacon, Moshinsky and Sharp and of Kemmer, 
Pursey and Williams are very often considered as 'classics' on the group theory of the collective model. Recently 
Rowc, Turner and Rcpka [49] gave an useful algorithm for the computation of SO (5) spherical harmonics. 
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A. Rowe and Bahri's work 

An alternative solution of the harmonic oscillator and of the Davidson potential has been given in ref. [17]. It is easy 
to recognize that the generalized coordinates of the nuclear collective model, q u and their conjugate momenta tt v , may 
be used to form operators that are closed under commutation. Having defined the scalars products 1 = ^„ | q y | 2 
and 7r 2 = J2u | 2 wc can introduce the operators 

Zi = 7T 2 Z 2 = (i 2 Z 3 = i £)(g„ • 7r„ + 7r„ • g„) (152) 

that span an sp(l,R) ~ 1) algebra. They are in fact closed under commutation: 

[Z u Z 2 ] = ~4iZ 3 [Z 3 , Z 2 ] = -2iZ 2 [Z 3 , Z x \ = 1iZ x . (153) 

With the linear transformation 

Xi = ±(Zi-Z 2 ) X 2 = iz 3 X 3 = - A {Z X + Z 2 ) (154) 
one may recognize the standard su(l, 1) ~ so(2, 1) commutation relations 

[X u X 2 } = -iX 3 [X 2 ,X 3 ]=iX 1 [X 3 ,X 1 }=iX 2 . (155) 

It is also very useful to the define raising, lowering and weight operators for this algebra (the so-called eigenoperator 
decomposition) 

X±=X 1 ± iX 2 A = A 3 (156) 

that obey the following commutation relations 

[A+, A_] = -2A [A ,A±] = ±A± . (157) 

The action of the above operators on orthonormal bases states for the irreps of su(l, 1) (| n, A) with n=0,l,2,..) is 
given by the equations: 



The Casimir operator is 



A+ | nX) = y/(X + n)(n + l) \ n + 1, A) 

A_ | n + 1, A) = i/(A + n){n + l)\ nX) (158) 
A | nX) = -(X + 2n)\ nX) 



X 3 — X\ — X 2 — A (A — 1) — A + A_ 



l(z x Z 2 + Z 2 Z^)-?f (159) 



and takes the values 



C | nX) = |a(A - 2) | nX) . (160) 

From what we have summarized here it follows that the harmonic oscillator hamiltonian may be written as (where 
we are omitting some hu factor) 

Z x +Z 2 = 4A 3 (161) 
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which is diagonal in the basis given above. Its spectrum is 

E nX = (2n + A) (162) 
where X = v + 5/2. With the remarkable nonlinear transformation [17, 53] 

h 2 e 

Z X ^Z X + — Z^^Z^ Z 3 ^ Z 3 (163) 

Z 2 

the new operators satisfy again the sp(l,R) ~ su(l, 1) commutation relations. The spectrum 



E nv = (2n + l + y/(v + 3/2) 2 + e) (164) 

is found introducing the definition of A that comes from the comparison between the eigenvalues of C obtained from 
(160) and (159). 

B. Algebraic approach to Coulomb-like and Kratzer-like potentials 

The preceeding treatment of the Davidson potential, together with ch. 18 in ref. [44], has served as a source of 
inspiration for the treatment given in [33]. The spectrum of the Coulomb-like and Kratzer-like potentials may be 
derived in a similar fashion by noticing that the following operators 

Zi = 4/3^ + ^) Z 2 = (3 Z 3 = 2(f-#-i) (165) 

are closed under commutation with the same relations of (153). With the Kratzer-like potential (that contains also 
the Coulomb-like case, when B = 0) the operator (3TL is in fact expressible as a linear combination of the elements of 
the algebra of su(l, 1), namely in the form 

0H = ZJA-A. (166) 

By defining again raising, lowering and weight operators we can write the eigenvalue equation for the Bohr hamiltonian 

as 

'(1 + 4e)X 1 - 2A+ (1 -4e)X 3 }y = (167) 

and following the procedure in [50] we can perform a (1,3) hyperbolic rotation of an angle to diagonalize the 
eigenvalue equation. By choosing tgh(6) = —(1 + 4e)/(l — 4e) (valid for e < 0) we obtain a diagonal relation: 

i 3 f = -4=* i ( i6s ) 

\J— 4e 

where ^ is the rotated wavefunction. The Casimir operator of the so(2,l) algebra is evaluated to be: 



--2 



= A z + X- - X+ + B + 2 , (169) 



with eigenvalue t(t + 3) + B + 2. The two last equations must be compared with the two following eigenvalue equations 
(for unitary representations D + [51]): 

x 3 \^0 = (£-</>) 1 0,0 

C 2 | </>,£) - 0(0+l)|0,O. (170) 
This comparison yields a spectrum of the form: 

A 2 1 A 

erf = 1 -. (171) 

,£ (V(t + 3/2)2 + S + 1/2 + 2 

that coincides with the one found from the direct solution of the differential equation with a Kratzer-like potential. 
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FIG. 28: Rotational part of the spectrum of the rigid triaxial rotor (Davydov's model). 

The algebra associated with the SO(5) group plays the role of a degeneracy algebra [52], while the group SO(2,l) is 
associated with the spectrum generating algebra. For what concern either the problem considered here and the one 
in the previous subsection the chain of subalgebras that gives the labels of the set of orthonormal states {| ^raLM}} 
is given as [17, 19, 53]: 

51/(1,1) x 50(5) D 17(1) x 50(3) D 50(2) 
A r a £ L M 

where A is an SU(1,1) lowest weight and a indexes the SO(3) multiplicity. These basis diagonalize the above problems. 
We can thus state that the problem studied so far displays a SO(2,l)xSO(5) dynamical algebra. 




C. Quasidynamical SO(2) symmetry for triaxial nuclei 

In fig. 28 we plot the eigenvalue of the rigid triaxial rotor hamiltonian (see sect. V) as a function of 70. For 7 = 
and 7 = 7t/3 the projection of the third component of the angular momentum on the intrinsic axis 3 gives a good 
quantum number (K), while for 7 = 7r/6 the eigenvalue of the projection on the intrinsic axis 1 is a good quantum 
number (R). In the intermediate regions, none of them may be taken as a good quantum number. Moving from 
7 = 30° towards 7 = 0°, different groups of states may be classified into bands: a first band (0 + , 2 + , 4 + , ...) tend to 
the finite axial rotor values; a second band (2 + ,3 + ,4 + , ...) is identified by its behaviour when 7 — ► (in fig. 28 this 
group of states somewhat cluster around 7 ~ 10° — 12°); the beginning of a third band (4 + , 5 + , ....) is seen to escape 
to infinity at a quicker pace (exiting fig. 28 at around 7 ~ 20°). The experimental observation that a classification 
in (3 and 7 bands seems an almost universal feature of nuclear spectra reinforces this choice. The labeling with the 
K quantum number is often encountered in the literature, although for what we have said here it is not adequate. 
One may describe this situation in terms of a quasidynamical symmetry [54, 55] of a somewhat strange character: at 
7 = 0° the group SO(2) is a symmetry of the system, associated with K, while at 7 = 30° another SO(2) group is a 
symmetry of the system, associated with R, being the chain U(5)dSO(3) common to the whole sector < 7 < it/6. 
In the intermediate region < 7 < tt/6 the SO (2) symmetry is broken (badly broken in a 'classical' sense), but it 
must be noticed (see fig. 28) that the structure of the rotational spectrum present at 7 = 0° persists in the whole 
sector without being altered in a dramatic way. Only a smooth and slight change may be seen. On the other side 
the structure of the 'maximally' triaxial rotor at 7 = 30° persists also in the region around 7 ~ 20° — 30°. In the 
intermediate region these groups of states escape to infinity, as already said. It must be further noticed that the 
regions where the various states that comes from the axial rotor side are more affected is exactly the region where 
the states coming from the 7 = 30° triaxial rotor diverge. The strange character of this quasidynamical symmetry 
mentioned above is that (at variance with the case discussed by Rowe and collaborators [54, 55] , where a true phase 
transition was present between two exactly solvable limits associated with different symmetries and different group 
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structures) here we are dealing with a smooth transition between two limits which formally have the same underlying 
group structure, SO (2), and there is no critical point in between. Therefore we conclude that the use of a label that 
mimics the K quantum number, that retain the formal division in (3 and 7 bands typical of an axial rotor, is not 
only justified by the empirical observation that non-axial nuclei display the same classification in bands, but it is also 
supported in view of arguments based on a group theoretical approach. It is not clear at present if a quantization 
procedure around a tilted axis may help to further shed light on this aspect. 

It is understood, however, that K has been introduced as a convenient label (for basis states and eigenvalues) that 
comes into play in the branching rules from SO(5) to SO(3) [47, 49]. 

VII. RECENT DEVELOPMENTS OF THE COLLECTIVE MODEL 

We cannot completely neglect, although we will just touch the argument briefly, that the collective model has 
recently experienced new important developments that are connected with the various solutions that we have sum- 
marized in the present review or that furnish new strategics to get numerical solutions. These developments have in 
common the Streben to simplification and tractability. 

A. Caprio's simplified approach to the CM 

The more general hamiltonian of the geometric collective model contains a series expansion in terms in the surface 
deformation coordinates and their conjugate momenta. The Bohr hamiltonian corresponds to the truncation of all 
the kinetic terms with order higher than the second. Caprio studied a simplified version [56] of the hamiltonian in 
which only the leading order kinetic term and the three lowest leading order potential terms are retained, namely: 

ff = ^M tt](°) + ^/3 2 - ^§C 3 /? 3 cos3 7 + ^/3 4 . (173) 

This hamiltonian has a structure rich enough to encompass rotational, vibrational and 7— unstable cases, but it's 
still very complex and contains four parameters. Caprio exploited analytic scaling relations to reduce the number of 
parameters from four to two and summarized the predictions of the model on two-dimensional contour plots. He used 
two 'simple' observations to make the reduction. Firstly, an overall multiplication of the hamiltonian by a constant 
does not affect eigenvalues and wave functions. Secondly, with the transformation 

V{J3, 7) V'(P, 7) = a 2 V(a/3, 7) (174) 

all the eigenvalues are multiplied by a 2 and the wave functions are radially dilated by 1/a. Other simplifications and 
graphs are given in the text to which we refer the reader [56] for a more detailed discussion. The author mentions 
also the difficulty in numerical diagonalization of the problem in the basis of the five-dimensional harmonic oscillator, 
emphasizing that a large number of basis functions must be used in order to have correct results. A possible solution 
to this problem is outlined in the following section. 

B. Rowe's tractable version of the CM 

The Bohr-Mottelson collective model may in principle be solved numerically for any potential that is an analytic 
functions of the two variables, but the expansion of rotational wave functions in a spherical vibrational basis is very 
slowly convergent. The diagonalization is therefore very difficult. Rowe proposed a method in which the basis is 
constructed with analytic wave functions whose p— part are centered around a non-zero value [57]. The convergence 
with this method is much faster for a wide variety of model hamiltonians and this method represents an important 
improvement over older models. This approach, that we will only briefly sketch, is formulated in an algebraic way using 
the group SU(l,l)xSO(5) as a dynamical group and recognizing that SO(5)- invariant hamiltonians are diagonalizable 
using SU(1,1) as a spectrum generating algebra (see [57] for details). The dynamical subgroup chain employed is 
SU(l,l)xSO(5)c SO(5)c SO(3). The Davidson potential [14-17], that is an exactly solvable case, is employed for 
what concern the dependence on /3. Its exact wave functions are (approximately) centered around the minimum of 
the potential well and give a good starting point for diagonalization. This simple version may be applied directly to 
all the cases which lie between the spherical vibrator and the Wilets and Jean (7— flat) limit. 
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TABLE VI: Low-lying eigenstates of the model described in [58]. Their total angular momentum is also indicated. A few other 
cases might be solved. 

A more general extension of use of the basis obtained from the solution of the Bohr hamiltonian with the Davidson 
potential is also treated. A periodic potential in gamma V(j) — — x cos 37 is added to the hamiltonian. The reader 
is referred to the original paper for the details and the discussion. This model is very flexible because of the rich 
variety of spectra that may generate and it is more physical because of the relaxation of the hypothesis of rigidity. 
Furthermore it allows for a tractable diagonalization procedure. 

In principle this model, with suitable choice of potentials and parameters, should be able to cover the whole 
complexity of nuclear collective features associated with the quadrupole degree of freedom. 

C. Rigid triaxial model with independent inertia moments 

In the rigid triaxial rotor model irrotational flow moments of inertia have been usually assumed. In Ref. [58] the 
rigid triaxial rotor model is studied including the three components of the moments of inertia as parameters that may 
be fitted to experimental data. The rigid rotor hamiltonian is 

3 

H = ^^Q-, (175) 

i=l 

where the three parameters Ai are defined as 

A 1 = -^, i = l,2,3 (176) 

Z-Li 

and I\ > I2 > 3-3- The operators Qi are the components of the angular momentum in the body- fixed frame of 
reference. The hamiltonian (175) may be rewritten in terms of diagonal and non-diagonal operators (using the 
definitions Q± = Q\ ± 1Q2) as 

H = AQ 2 + FQ 3 + G(Q 2 + + Q 2 _) (177) 

where A = (Ai + A-i)f2, F = A3 — A and G = (A\ — A-i)f\. A few low-lying eigenstates may obtained analytically 
from the diagonalization of the hamiltonian in the rigid rotational basis | LMK) and they are summarized in Table 
VI. Some of the eigenvectors and electric quadrupole properties are also discussed in the paper and from an analysis 
it emerges that the components of the moment of inertia are not irrotational. For further discussion we refer the 
reader to the cited work. 
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APPENDIX A 



Number of repetitions 



Within a given set of quantum numbers (a given r in the general 7— unstable case, or N in the harmonic oscillator 
case) the same L may appear more than once. We define the repetitions of L as the number of counts of that 
L minus one. Despite the fact that this appendix could seem a purely academic divertissement, its novel number 
theoretical content is not completely detached from its physical significance: if a model predicts two states that are 
degenerate in energy to have the same L quantum number it may be very difficult to experimentally disentangle the 
two. The procedure of assignation of quantum numbers to measured spectra is not free from its own uncertainties 
and the outcome is a crucial test for theories. It seems thus very important to have a detailed knowledge of what we 
should expect from our theoretical calculation to look through the data. It will be clear that, in the following cases, 
the occurrence of repetitions is rather rare and it appears frequently only for high lying states so that likely is not 
prejudicial to any analysis that was carried on so far. 



General 7— unstable case 



We will begin with the general 7— unstable case: the problem is to count how many repeated multiplicity of SO (3) 
appears within a given irrep of SO (5), or to count the repetitions of L's for a given r. Looking at table I in the text 
we can see that up to r = 5 every L appears just one time, so the number of repetitions is 0. For r = 6, the quantum 
number L = 6 appears two times, so the number of repetitions is 1. Continuing along this line one easily finds (by 
hand in the first few cases) that the number of repetitions, starting from r = 6 is: 



1, 1, 2, 4, 5, 7, 10, 12, 15, 19, 22, 26, 31, 35, 40, 



(A.l) 



This sequence has been calculated with a simple computer code and it corresponds to the sequence obtained from 



R(r) 




(A. 2) 



with n = t — 3 = 3, 4, 5, where square brackets denote the integer part, or floor function. Within the same range 
of n this is also connected to the number of solutions to the equation x + y + z = ( mod n ) with 0<x<y<z<n. 
For instance if n = 3 the only possible solution is x = 0, y — 1 and z = 2. This sequence arises also in connection 
with the so-called 'orchard problem' (See [59]). 



Harmonic oscillator case 



For the case of harmonic oscillator a further analysis may be done. For a given iV a certain number of pairs of 
(nfi-, t) satisfy N = 2np + r. For each pair a value of r is associated with a set of quantum numbers L (the same of 
the preceding section). Hence the set of L's corresponding to each ./V is larger that in the former case (the degree of 
degeneration is higher) and the sequence of numbers of repetitions, starting from N = 4, may be written as follows 

2, 3, 7, 11, 18, 26, 36, 48, 63, 79, 99, 121, 146, ... (A . 3) 

while in all the cases with TV < 3 there are no repetitions as can be seen from Table II. Unlike the previous case it 
was not possible to make a connection between the present sequence and any previously known sequence [60], nor it 
was possible to find a simple formula to express the Ath term of the sequence. Note that the number of repetitions 
in this case is not simply given by 

[N/2] 

R(N) = ^ R(t = N- 2np) (A . 4) 

because it happens that many repetitions come from different (n/3,r) pairs (otherwise they would have started from 
N = 6). It is furthermore intriguing to notice that when the terms of the sequence (A . 3) are subtracted with the 
corresponding terms of the sequence that may be built from formula (A . 4), the remaining sequence (starting from 
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N = 4) is 

2, 3, 6, 10, 15, 21, 28, 36, 45, 55, 66, .... (A . 5) 

and apart from the first two elements the differences between two of any other couple of adjacent terms increases 
arithmetically (!). 

Two very simple computer codes, IRREP and IRHAR, written in Fortran77, are available from the author or from 
WWW [61]. They count the number of times that a given quantum number occurs, and give thus the number of 
repetitions as in (A . 1) and (A . 3) . 

APPENDIX B 
Meaning of the terms "softness" and "rigidity" 

The purpose of the present Appendix is to clarify the meaning that should be attributed to the words "softness" 
and "rigidity" in the context of the collective model. 

In general the solution of the Schrodingcr equation depends on the potential that is chosen. In the case of the Bohr 
hamiltonian (restricted to the quadrupole degree of freedom) the potential is a function of < (3 and < 7 < 7r/3. 
Every couple of values (3, 7 univocally corresponds to a given quadrupole surface. The state of the quantum system is 
represented by wavefunctions that have a dependence on (3, 7 or both. When the potential is very deep and narrow its 
wavefunctions (especially the more low-lying ones) are very localized in the region of the minimum and they may be 
associated with just one value of (3. The object (that lives in the configuration space) that is described by such state 
has therefore a rigid shape in the sense that its surface doesn't exhibit fluctuations. On the contrary a more shallow 
potential gives rise to extended wavefunctions. This means that the real object may be thought as a superposition of 
different quadrupole surfaces each with its own weight. The term "soft" was used to describe such an object. This 
system exhibits non-null fluctuations around the mean square deformation. 

Sometimes it is argued by analogy with the work of Wilcts and Jean [11] that the independence of the potential on 
one variable is associated with the term softness. The dependence or independence on a given variable should not be 
used to establish the softness or rigidity of some state. The unstable case of Wilets and Jean (for what concern the 
7 d.o.f., that is a sort of angular variable and not a radial variable) is, among the many possible soft cases, the most 
extreme one. 

For example wc may decide to set the potential in (3 to a constant on its whole semi-infinite domain. 

In this f3— independent case the cigenstatcs lie in the continuum and the eigenfunctions have a periodic behaviour at 
infinity. It is thus not possible to associate a finite value to the deformation of such a state! The situation in the paper 
of Wilets and Jean is different: two cases were discussed, the infinite square well and the displaced harmonic oscillator. 
They are both examples of a soft behaviour, because the two quantum mechanical potentials give rise to solutions 
endowed with extended wavefunctions. One can see that while the displaced oscillator has a true dependence on /3, 
the infinite square well has a dependence hidden in the boundary conditions (that are essential to fix the spectrum). 
The analogy with the 7— unstable case is thus rather feeble, because a f3— unstable counterpart has not a definite 
physical meaning. The word "unstable" refers to the fact that the potential has no minima (and hence there's no 
"stability" in the classical sense) and should be read as a synonymous of "independent". 

In summary the word "soft" applies to the usual quantum mechanical case, while the word "rigid" must be associated 
with very peculiar cases for which either the value of (3 (and/or 7) is fixed or the potential is of a very special kind. 
At the beginning of the study of the Bohr-Mottelson model many "rigid" cases were discussed, mainly for the need 
of simplicity and for the lack of more refined ("soft") solutions. 

APPENDIX C 

Electromagnetic transitions 

In some of the spectra displayed in the paper, electromagnetic transition rates have appeared. Although it is not the 
subject of the present monograph, it is useful to discuss to some extent how they may be derived, especially because 
BE(2) values, being very sensitive to the details of the wave functions, are often considered as crucial quantities for 
the application of the aforementioned models to nuclei. 
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The electromagnetic operator that is often encountered in this context reads: 

T{E2) = t/3 



D$ cos 7 + 7^(^ 2 2 + <- 2 ) sin 7 



(C.l) 



where t is a scale factor. This is the first order approximation of the more general E2 operator [4, 63]. In the geometric 
collective model this operator is written in the laboratory frame coordinates as an expansion: 



T(E2; n) = aa 2 , u + b[a x a]< 2) + 0(a 3 ) , 



(C.2) 



where each term is coupled to a total angular momentum equal to 2. It has been often argued [32, 62] that a second 
order approximation 



T(E2) = af3\D^ cos 7 + -j=iP™ + D%_ 2 ) sin 7 



(2) 



-,(2) 



cos 2 7 --^(D M)2 



1 



(2) _ LjD (2^) sin 2 7 



(C.3) 



leads to non-trivial modifications in the calculations of the matrix elements. Transition rates are then calculated in 
the following way: 



B{E\; Ji — ► Jf 



(J, | T{E\) | J f ) 
2Ji + l 



(C.4) 



where A is the angular momentum that characterizes the transition. Usually B(E2) values have been normalized 
taking the value of the B(E2) for the transition 2^ — > Of as 100. 

Finally we mention that, for E0 transition, the electromagnetic operator is 



T(E0) oc (3 2 



(C.5) 
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